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ABSTRACT 

We use recent observations of the He I /U0830A absorption line and 3D hydrodynamical 
numerical simulations of the winds collision, to strengthen the case for an orientation of the 
semimajor axis of the massive binary system r] Carinae where the secondary star is toward us 
at periastron passage. Those observations show that the fast blue absorption component exists 
for only several weeks prior to the periastron passage. We show that the transient nature of 
the fast blue absorption component supports a geometry where the fast secondary wind, both 
pre and post-shock material, passes in front of the primary star near periastron passage. 

Key words: Stars: individual: rj Car - Stars: winds, outflows - Stars: mass-loss - (Stars:) 
binaries: general - Infrared: stars 



1 INTRODUCTION 

T] Car is a very massive stellar binary system, with an orbital pe- 
riod of 5.54 yr (Damineli 1996), as observed in all wavelengths 
(e.g., radio, Duncan & White 2003; IR, Whitelock et al. 2004; vis- 
ible, van Genderen et al. 2006, Fernandez-Lajus et al. 2009; UV, 
Smith et al. 2004; emission and absorption lines, Nielsen et al. 
2009, Damineli et al. 2008a, b; X-ray, Corcoran 2005, 2010, Ham- 
aguchi et al. 2007). The high eccentricity of e ^ 0.9 results in rapid 
changes in emission and absorption lines, as well as in the contin- 
uum, near each periastron passage. The several weeks of the rapid 
changes occurring every orbital period is termed the spectroscopic 
event. These lines might originate in different places in the binary 
system: the primary star, the secondary star, their respective winds, 
and the colliding winds structure which is termed the conical shell. 

As tj Car is the best studied binary luminous blue variable 
(LBV), it holds the key to our understanding of other LBVs. It 
is particulary important to understand the behavior near periastron 
passage, where the strongest binary interaction takes place, and for 
that is crucial to know the orientation of the binary system. Namely, 
the direction of the primary more massive LBV star relative to its 
less massive but hotter companion at periastron passage. While it is 
agreed that the inclination of the binary system is i — 41° (David- 
son et al. 2001; Smith 2006), there is no agreement on the direction 
of the semimajor axis, termed periastron longitude. The orientation 
is measured by the angle a>: a> = 0° for the case when the secondary 
moves toward us before periastron passage and the semimajor axis 
is perpendicular to the line of sight, to = 90° when the secondary 
is toward us at periastron passage, and to = 270° when the primary 
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is toward us at periastron passage. Several properties of the binary 
system have been used to deduce the orbital orientation, with con- 
tradicting results (for details see Kashi & Soker, 2008b, 2009c). 

One of the properties that led to a contradicting conclusion on 
the orientation is the behavior of the blue absorption wing of the 
He I il0830A line. In Kashi & Soker (2009b) we constructed a 
toy model where the material responsible for the blue absorption 
wing was assumed to reside in the colliding winds region - the 
conical shell - close to the binary system. This model is able to ac- 
count for the transient appearance of the blue absorbing wing and 
to the finding that the maximum absorbing velocity is reached sev- 
eral days before periastron passage, only if the secondary is toward 
us near periastron passage, i.e., to ^ 90°. In a recent paper Groh et 
al. (2010; hereafter G2010) reached an opposite conclusion. Com- 
paring their observations with a model based on 3D numerical sim- 
ulations of the colliding winds structure, G2010 suggested that the 
orientation is that of to = 243°. Namely, the primary is closer to us 
just before periastron passage. In this paper we critically reexamine 
both models. 



2 PROPERTIES OF THE ABSORPTION PROFILE 

The He I /U0830A line has been observed for more than a decade, 
but only occasionally (Damineli et al. 1998, Groh et al. 2007; 
Damineli et al. 2008b; G2010). The line shows a P-Cygni profile 
varying in time, especially close to the spectroscopic event. It has 
a complicated emission with three spikes, and absorbing compo- 
nents with velocities of up to -v - 1900 km s" 1 . In this section 
we describe the relevant properties of the blue absorption wing that 
rapidly appears and disappears near periastron passages. 

G2010 observed the He I il0830A across the 2009 periastron 
passage. In Fig.[T]we show the absorption profile at phase 11.998, 
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Figure 1. The flat part of each of the three absorption components identified 
in the paper , marked on the absorption profile at phase 11.998 (four days 
before the periastron passage of the 2009 spectroscopic event) taken from 
G2010. 



four days before the 2009 spectroscopic event, assumed to occur 
at phase 12.000 (Damineli et al. 2008b, Groh & Damineli 2004), 
taken from fig. 1 of G2010. The profile at phase 11.991 is practi- 
cally similar in shape, but the absorption is somewhat weaker. 

At phases 11.998 and 11.991 we can identify three flat parts 
in the spectrum, as marked on Fig.Q] 

(i) P-component. The deepest component can be safely iden- 
tified with the regular primary wind. This part exists during all 
phases, it is deep, and its velocity matches that of the primary wind. 
The primary wind reaches terminal velocities of 500-650 km s -1 , 
depending on latitude (Smith et al. 2003). There is no dispute about 
the origin of this component. 

(ii) M-component. Although less prominent, the M-component 
has a flat profile in the velocity range -v = 620-720 km s -1 . We are 
not sure about its nature, but we do see signatures that it exists at 
phases 11.875, 12.014 and 12.041, although less prominently. We 
attribute it to parts in the primary wind that are accelerated to some- 
what higher velocities. These may be, for example, clouds with 
higher opacity, e.g., due to dust formation. We note that before the 
2009 periastron passage the primary wind might have undergone 
some fluctuations in its velocity, which may have been evident in 
the X-ray lightcurve (Kashi & Soker 2009d). 

(iii) F-component. We identify the main parts of this component 
with the flat part in the velocity range -v ^ 1150-1400 km s -1 , 
and the tail formed by faster and lower column density gas with 
absorption extending to -v = 1900 km s _1 . This tail is not due 
to line broadening in the main part of the absorbing component. 
The F-component does not exist at all at phases 11.875, 12.014 
and 12.041. Namely, it is a transient component appearing for a 
short time before periastron passage. In its flat part this component 
absorbs ~ 25 - 30 per cent of the continuum at phase 1 1.998, and 
~ 20 per cent at phase 1 1.991. 

We attribute the F-component to pre and post-shock secondary 
wind material close to the binary system (Kashi & Soker 2009b). 
The post-shock gas resides in the conical shell. We further assume 
that the secondary wind gas absorbs most of the continuum emitted 
by the gas in the conical shell (emission that increases just before 
periastron passage), in addition to absorbing some fraction of the 
continuum emitted by the primary wind. It is this component that is 
in dispute between G2010 and us. Although G2010 and us attribute 



it to the colliding winds, in the model of G2010 the absorbing gas 
is at a distance of tens of AU from the binary systems, while in 
our model the absorbing material is relatively close to the binary 
system. 

One thing to notice is that the first observation after the 2009 
spectroscopic event is at phase 12.014, before the system exits from 
the X-ray minimum (Corcoran 2010), and before the hard ionizing 
radiation from the secondary resumes after exiting the dense pri- 
mary wind. As we note in section [3] this holds difficulties to the 
model of G20 10. 



3 PROBLEMS WITH THE MODEL OF G2010 

The purpose of this section is to set the motivation for a model 
where the fast absorbing material reside close to the secondary star, 
and the secondary star is in front of the primary near periastron 
passage. We critically examine a model with an opposite binary 
orientation, and reveal its problems. 

G2010 present (in their fig. 12m-o) the total column density 
as a function of gas velocity at 5 phases according to their assump- 
tions. Their main assumptions are that: (1) The photospheric radius 
of the primary at 1.083 yum is 2.2 AU. (2) A point source can be 
used for the central source of the 1.083 [im continuum emission. 
(3) The central source is the only source of the continuum (below 
we claim that the conical shell is an additional source). (4) G2010 
did not consider even qualitatively the change in level population. 
In their model the change in level population must be significant; 
below we argue why this is unlikely. 

In more details, their fig. 12 gives the density of material as 
a function of line of sight distance to the primary for 5 orbital 
phases (11.875, 11.991, 11.998, 12.014 and 12.041) and 6 orien- 
tations (en = 0°, 50°, 90°, 180°, 243°, and 270°). In addition, the 
column density along the same assumed line of sight to the pri- 
mary star is given as a function of the line of sight velocity, for 
each orientation and phase. Fig. 12 of G2010 allows us to deduce 
the basic properties of their model. 

G2010 did not relate the total column density to the absorp- 
tion of the He I /110830A line. We attempt to do so and check what 
the results of G2010 imply. As we now show, their model, with a 
preferred orientation of a> = 243°, suffers from severe problems 
because the column density does not follow even qualitatively the 
behavior of the observed absorption. In their attempt to resolve this 
discrepancy, G2010 claim that the presence of high velocity gas 
in our line of sight is a necessary, but not sufficient condition for 
the presence of high velocity absorption. As the amount of absorp- 
tion depends on the population of the lower energy level of the 
He I /110830A line (2 3 5), they attribute the differences in absorp- 
tion to different population of the 2 3 5 atomic level. 

3.1 Pre-periastron passage 

For the three phases before periastron passage (11.875, 11.991, 
11.998), G2010 model (in their fig. 12o) predicts the absorbing 
column density in the range -v = 2000-3000 km s to be about 
equal to that in the velocity range — v = 1600-1800 km s . Obser- 
vations show substantial absorption in the -v = 1600-1800 km s _1 
range at phases 11.991 (~ 3 - 10 per cent absorption) and 11.998 
(~ 5-13 per cent absorption), while there is no absorption at 
-v > 1900 km s 1 (fig. 1 of G2010). G2010 attribute the differ- 
ence to stratification in the population of the 2 3 S level, from which 
the He I /110830A line is absorbed. However, they do not present 
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a qualitative study of the way the population of the 2 3 S level can 
evolve so rapidly. We find it extremely unlikely (unless there is a 
carefully arranged fine tuning) that differences in populating the 
2 3 S level can account for this difference. The reason is as follows. 

According to G2010 model during the phase interval 1 1.991- 
11.998 the absorbing gas at — v = 2000 km s _I resides at a dis- 
tance of r a - 3.5a from the binary system, where a - 16 AU is 
the semimajor axis (their fig. 12n). They also find that the gas ab- 
sorbing in the velocity range -v = 1600-1800 km s -1 resides at 
r a =i 2.5 - 3a. During this phase interval, the binary orbital sep- 
aration is r < 0.2a. Therefore, the differences in the radiation of 
the two stars at r - 2.5 - 3a and r — 3.5a are small, and can- 
not account for the required large differences in the population of 
the 2 3 S level. Moreover, the increase in the amount of absorption 
as the system approaches periastron in the model of G2010 is at- 
tributed to an increase in the population of the 2 3 S level. This is 
probably (G2010 do not specify the reason) because the secondary 
gets deeper into the primary wind, and less of its radiation reaches 
the absorbing gas. It is expected that the outer regions, where gas 
moving at -v = 2000 km s resides, will get less of the secondary 
radiation, and hence will absorb more. This is opposite to what their 
model requires. 

A potential way to handle the problem is that the material 
closer-in is of higher density, and it recombines first or its colli- 
sional excitation is more efficient. This would not work either. The 
rate of these two processes is linear with density. The density ratio 
between the region absorbing in the -v = 1600-1800 km s -1 range 
to that absorbing in the -v = 2000 km s range is ~ 2 (G2010; fig. 
12m,n). Therefore, the absorption at -v - 2000 km s -1 is expected 
to be one half of that in the -v = 1600-1800 km s -1 velocity range. 
This amounts to ~ 1.5-5 per cent of the flux at phase 1 1.991 and 
~ 2.5 - 6.5 per cent of the flux at phase 11.998. This is clearly 
ruled out by observations, as the absorption at -v = 2000 km s -1 
observed in these two phases is practically zero. 

3.2 Post-periastron passage 

In phase 12.014 (28 days after phase 12) the column density calcu- 
lated by G2010 in the velocity range -v = 1200-1400 km s" 1 is a 
factor of 2-3 lower than that at phase 1 1.998, while in -v = 1600- 
1800 km s _I this ratio is < 2 (their fig. 12o). At phase 12.014 the 
secondary radiation has not started yet to ionize the material at large 
distances, as it is still inside the dense region of the primary wind. 
We know this from the behavior of the highly ionized lines, e.g., 
in the Weigelt blobs, which show that the hard radiation resumed 
(in the 2003.5 spectroscopic event) well after the X-ray minimum 
ended (Damineli et al. 2008b; Mehner et al. 2010). In the 2009 
spectroscopic event the X-ray minimum ended at phase ~ 12.014 
(Corcoran 2010), when one of the He I /110830A line measure- 
ments was conducted. Namely, at this phase the secondary hard 
ionizing radiation did not ionize the He in the relevant distance of 
tens of AU. The density in the relevant region at phase 12.014 is 
not much different from that at phase 11.998. If it was for recom- 
bination, there should have been more time for the gas at phase 
12.014 to recombine. 

Over all in the velocity ranges —v = 1200-1400 km s and 
-v = 1600-1800 km s _1 the column density calculated by G2010 
at phase 12.014 is ~ 0.3 - 1 that at phase 1 1.998, and the ionization 
state should be the same. Therefore, it is expected that the absorp- 
tion at phase 12.014 in the velocity ranges -v = 1200-1400 km s -1 
and -v = 1600-1800 km s _1 would be -0.3-1 times that at phase 
1 1.998, or even above its value at phase 1 1.998. This is contrary to 



observations, as there is no absorption at all for -v > 1 100 km s 
at phase 12.014 (fig. 1 of G2010). This is a severe problem for the 
model presented by G2010. 

Not only the model of G2010, but any other model that would 
assume absorption by materiel residing at a large distance (a few 
xlOa or more) is problematic. The reason is that the absorption of 
the He I /110830A line in high velocities is a transient event, which 
occurs very close to periastron, and must be related to variations 
in the system shortly before periastron. The material far from the 
binary system does not suppose to show fast variations close to 
periastron, and for that it is a bad candidate for an absorber of the 
HeI/U0830Aline. 



4 THE MODEL 

4.1 The column density of the absorbing gas 

We perform hydrodynamic simulations of the colliding winds to 
show that the column density of the conical shell and the pre- 
shocked secondary wind at high velocities is high enough to ac- 
count for the absorption of the He I /110830A line. The simulations 
are performed with Virginia Hydrodynamics-I (VH-1), a high reso- 
lution multidimensional astrophysical hydrodynamics code devel- 
oped by J. Blondin and collaborators (Blondin et al. 1990; Stevens 
et al., 1992; Blondin 1994). The code includes radiative cooling at 
all temperatures of T > 2 x 10 4 K. The radiative cooling is set to 
zero for temperatures of T < 2 x 10 4 K, and the gas cannot cool 
to lower temperatures. The cooling function A(T) (for solar abun- 
dances) is taken from Sutherland & Dopita (1993; their table 6). 

The numerical simulations were performed in the Cartesian 
geometry (x,y, z) mode of the code (a 3D calculation), where the 
orbital plane is the ( x, y) plane with the x axis taken to be parallel 
to the line connecting the two stars. Our numerical simulation grid 
has 112 equal-size grid points along each axis. The primary and the 
secondary are located in the middle of the y and z axes. The binary 
separation between the two stars is set to be constant at each numer- 
ical run (for the justification see Akashi & Soker 2010). The mass 
loss rates and velocities of the winds are Mi = 3xlO" 4 M G yr' and 
M 2 = 10~ 5 M Q yr -1 , and vi = 500 km s _1 and v 2 = 3000 km s~\ 
respectively (e.g., Pittard & Corcoran 2002; Kashi & Soker 2008b 
and references therein). We use eccentricity of e = 0.9 and semi- 
major axis of a = 16.64 AU. We performed four runs with binary 
separation of r = 2, 3, 4, and 10 AU, and in each created a density 
profile and wind velocity field around the two stars. At the location 
of each star we inject its appropriate wind. We let the flow reach an 
approximately steady state. An example of the results of the simu- 
lation is shown in Fig. [2] 

For each run we calculate the direction to the observer for a 
system with (co, i) = (90°, 41°), and for the phase before periastron 
corresponding to the binary separation of the run. We take a plane 
perpendicular to the line of sight and located behind the primary, 
and from this plane we calculate the column density along a grid 
of lines of sight to the observer. We select the location of this plane 
in an approximately optimal way to include as much of the grid as 
possible. This procedure favors accuracy in the calculation of the 
column density of the high velocity gas towards the observer. Each 
line of sight is divided into cells, as the number of z planes the line 
of sight crosses in our grid. We calculate the column density along 
each line of sight as function of the velocity projected along the 
line of sight. Like G2010 we use velocity bins of 50 km s . 

We adopt the conclusions of G2010 that the main source of 
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Figure 2. Our hydrodynamical simulations showing density and velocity 
maps for an orbital separation of r = 2 AU. The bar on the right gives the 
density color-code in units of log[p (g cm" 3 )]. The thin arrows represent the 
velocity in the (x,y) plane (the orbital plane) and are scaled with the axes 
such that one unit represents 5000 km s . The secondary is the star on the 
right at (x, y) = (5, 4) AU. The arrows are in the line of sight direction to 
an observer at (a), i) = (90°, 41°) and at phase 0.9964. The column density 
from each line of sight is weighted with a Gaussian (equation[T}. 



the 1.083 fim continuum emission comes mostly from a region of 
size ~ 2.2 AU centered on the primary star. We term it the central 
source. Though G2010 used a point source for the central source 
of the 1.083 yum continuum emission, they mentioned that it is of 
size of several AU, based on observations by Weigelt et al. (2007) 
in the K-band, which were calibrated to 1.083 yum (which is within 
the I-band), with the help of the theoretical results of Hillier et al. 
(2001). Weigelt et al. (2007) found the decrease function of the cen- 
tral source intensity with projected distance from the center r p to be 
close to a Gaussian (their fig. 6). We therefore take the emission to 
be from a radius of R cm = 4.4 AU, and weighted the intensity as a 
Gaussian with a = 2.2 AU 

I = I Q exp(-r^/2cT 2 ); for r p < R era , (1) 

where I is the intensity at the center (the cut-off at R em is used due 
to numerical constrains and results in a negligible error). 

The column density from each line of sight is weighted with 
the Gaussian in equation 10 according to its distance r p from the 
line of sight of the primary (Fig. [2}. This procedure gives us an 
effective column density in each velocity interval and at each phase. 
The results are presented in Fig. [3] 

From Fig. [3] we learn the following: (1) For a> = 90° the 
column density in the conical shell and the undisturbed wind is 
high enough to account for the absorption of the He I /U0830A 
line for our favored orientation of 0) = 90°. As the model of 
G2010, our model is unable to calculate the population of the 2 3 S 
level of He I, and hence we cannot calculate the exact absorp- 
tion in the He I /110830A line. (2) The high column density for 
-v > 1000 km s -1 is a transient component which appears only for 
a short period before periastron passage. (3) The column density 
clearly increases as the system approaches periastron (note that in 
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Figure 3. The calculated column density as a function of projected velocity 
along the line of sight. 



fig. 12 of G2010 this is not the case). (4) The column density de- 
creases towards higher blue velocities, despite of minor fluctuations 
due to the limited numerical resolution. (5) At short binary separa- 
tions (r < 4 AU; t > 0.987) there is a plateau in the column density 
between ~ 1400 km s -1 and ~ 2100 km s -1 , followed by a decrease 
in column density towards higher blue velocities. 

Our model gives that at very high velocities of -v > 
1900 km s the column density is also large, namely we have more 
absorbing gas than required to account for the observations. This 
discrepancy may be the result of an anisotropic secondary wind 
velocity close to periastron passage, that can occur if the accelera- 
tion of the secondary wind is not isotropic, an effect which was not 
taken into account in our simulation. Such an effect can reduce the 
amount of gas which outflows at -v > 1900 km s -1 , The accelera- 
tion can be efficient along polar directions, but below i ^ 60° radi- 
ation from the primary and previous accretion of blobs might cause 
the secondary wind to be slower. This explanation cannot work for 
correcting the surplus absorbing gas in the model of G2010, as their 
absorbing material is much further away from the colliding winds 
region. 



4.2 Covering the emission source with the conical shell 

From Fernandez-Lajus et al. (2010) we find that the I-band con- 
tinuum has a quiescence value of / = 3.55 mag. We sample the 
increase in the I-band continuum flux above the quiescence value, 
f c (t), from the observations of Fernandez-Lajus et al. (2010). These 
observations show that close to the 2009 periastron passage the I- 
band continuum flux (and therefore the 1.083 fim line flux) has in- 
creased by f c (t = 0) - 20 per cent above the quiescence value (Fig. 

We attribute most of this extra emission to the conical shell close 
to the binary system. We further assume that there is a strong ab- 
sorption by the conical shell and the pre-shocked secondary wind 
of the HeI/110830Aline. 

We analytically calculate the fraction of the central source of 
the 1 .083 /im continuum emission region that is covered by the con- 
ical shell (namely, covered behind the conical shell from the line of 
sight of the observer). The parameters of the stars and the conical 
shell are taken as in our previous papers (Kashi & Soker 2009b, c, 
and references therein). The contact discontinuity is approximated 
as an hyperboloid with semi-opening angle <p a si 60°. This is an 
approximation, as numerical simulations show that the structure is 
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Figure 4. The geometry of the binary system and the 1 .083 /im continuum 
emission region. 

much more complicated (Akashi & Soker 2010; Parkin et al. 201 1). 
We take into account the rotation of the conical shell close to peri- 
astron, caused by the non negligible orbital velocity relative to the 
primary wind velocity, and the acceleration of the primary wind 
(Kashi & Soker 2009b, c). The eccentricity is taken to be e = 0.9, 
such that the orbital separation at periastron is a,, = 1.66 AU. This 
implies that the secondary star and part of the conical shell are in- 
side the central 1.083 yum continuum emission region near perias- 
tron passage. Our assumptions break down after periastron passage, 
when the conical shell collapses. 

We again use our preferred orientation (co, i) = (90°, 41°) The 
geometry of the calculation is shown in Fig. [4] For every orbital 
phase we divide the observed half sphere around the primary that 
is closer to the observer into spherical volume elements (bins). For 
each emission source element we check whether the conical shell is 
between the observer and the emission source element, or whether 
it is within the conical shell itself. If the conical shell is between the 
observer and the emission source element, we use our assumption 
of high optical depth in the conical shell, and assume the emis- 
sion from the covered source is totally absorbed. In other words, if 
the emission source element is either of the two (behind the con- 
ical shell or inside it), we consider this element to contribute to 
the 1.083 ytim line absorption. The contribution of each element is 
weighed as the Gaussian, with the projected distance to the line of 
sight, as discussed above. 

The calculated flux of the 1.083 jim flux line according to our 
model is 

FobsW I - f„(t) + fc(t)d ~ fs) 
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Figure 5. Behavior of several properties just before periastron passage. The 
I-band continuum fluxes before the 2009.0 and 2003.5 periastron passages 
are very similar. Circles and squares :The observed decrease in the flat part 
of the F-component prior to the 2009.0 and 2003.5 periastron passages, 
respectively. The error bars are calculated from the maximum and minimum 
absorption values in the velocity range of the F-component. 



(the primary star). When to = 90° our model accounts for the fast 
increase in absorption of the flat part (F-component; Fig.QJ. Con- 
sidering the many uncertainties and the simplicity of our model, our 
results are in reasonable agreement with the observed decrease in 
the flat part of the F-component in the velocity range -v ^ 1 150— 
1400 km s _I (the tail to — v = 1900 km s -1 is part of this component, 
and behaves in the same way). 

The sizes of the central and extra emission sources are impor- 
tant for this modelling. G2010 use a point-like continuum source in 
their quantitative calculation of the column density of the absorbing 
material (total column density, not only that of the relevant atoms). 
For the model of G2010, where co = 243°, it does not make much 
difference to use a point source, as their absorbing material is sev- 
eral tens of AU from the source. However, calculating the column 
density to a point source fails to reproduce results of a model where 
the absorption occurs a few AU from the stars, such as our model 
(Kashi & Soker 2009b), which we improve here. G2010 do discuss 
(in their section 5.3) a 2.2 AU extended source for an orientation 
of 10 = 90°, and find the duration of the high velocity absorption 
to be too short. Our point here is that an extra source must be in- 
cluded to take into account the increase in the I-band and for the 
extra absorption. In our model the extra source is the conical shell. 
Our results show that one must use an extended emission source, 
both for the central source (the primary star) and the extra source 
from the conical shell. Using the two extended sources overcomes 
the objection of G2010 to the co = 90° orientation. 



et al. 2010, as defined above), f p is the fraction of the central contin- 
uum source covered by the conical shell, f c is the increase in the I- 
band continuum flux above the quiescence value (Fernandez-Lajus 
et al. 2010), and /, is the fraction of the conical shell emission that 
is absorbed by the F-component. If the conical shell self absorb all 
its emission, then /, = 1, while if it absorbs only from its volume 
that is behind the conical shell,then /, = 0.5. We present results 
for both these values that bound our expectation. Fig. [5] presents 
fobsW/fcontinuum(0 close to periastron, as well as f c and /,,. The 
conical shell absorbs up to f p (t = 0) = 25 per cent of the inten- 
sity of the 1.083 //m continuum emission from the central source 



5 SUMMARY AND DISCUSSION 

The He I /110830A line of 77 Car has been observed across the 2009 
periastron passage by G2010. We identify three main absorbing 
components in the He I /110830A line profile (Fig. [TJ. Two com- 
ponents can be identified with the dense primary stellar wind; the 
P-component (500-650 km s -1 ) and the M-component (— v = 620- 
720 km s -1 ). The third, more interesting F-component shows a flat 
part in the velocity range -v - 1150-1400 km s , with a tail to 
-v = 1900 km s . The tail results from faster moving gas, and 
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not from line broadening. There is no absorption in this line at 
-v > 1900 km s" 1 . 

G2010 present their interpretations of their observations, con- 
cluding that the orientation of the binary system is such that a> = 
243° (the primary closer to the observer just before periastron). In 
section [3] we point out some problems in the model of G2010, in 
particular that they cannot account for the transient nature of the 
line. The model of G2010 fails to reproduce, even qualitatively, the 
observed absorption profiles of the He I /110830A line, and can- 
not account for the absence of the F-component after periastron 
passage. In analyzing our (Kashi & Soker 2008b, 2009b) preferred 
orientation of a> = 90°, G2010 use a point-like central continuum 
source for calculating the column density of the absorbing material. 
Based on that, they claim that the <±> = 90° orientation is not com- 
patible with their observations. In section [4] we find that the usage 
of a central point source does not do justice to the a) = 90° case. 
This is the main reason why G2010 fail to reproduce results of a 
model where the absorption occurs a few AU from the stars, like 
in our model. 

We improve our pervious model (Kashi & Soker 2009b) which 
assumed that the conical shell is the main absorber of the high ve- 
locity component of the He I /U0830A line (section[4]l. Our model 
and results can be summarized as follows: 

(i) We assume that the central star is a Gaussian weighted ex- 
tended continuum source as given by equation |QJ. This assumption 
stands on a solid ground (G2010; Weigelt et al. 2007). 

(ii) We assume that the increase in the 1 .083 yum continuum near 
periastron passage comes from the conical shell (the collision re- 
gion of the two winds). The main source is the shocked primary 
wind. This is based on analysis we performed in previous papers. 

(iii) We assume that the main absorbing gas of the He I /110830A 
line is both pre and post-shock material in the fast secondary wind. 
This assumption stands on a solid ground. Based on recent numer- 
ical simulations (Akashi & Soker 2010), we showed in section [4] 
that the column density of the fast moving gas is as high as found 
by G2010. Some segments of the wind are at a temperature of 
~ 10 4 K, and there are enough recombining He atoms. This will 
be calculated in a future paper based on high resolution numerical 
simulations. 

(iv) We concentrate on the absorption of the F-component. We 
assume that in the flat velocity range of the F-component (see Fig. 
[T}, the absorbing gas in the conical shell absorbs most of the conical 
shell emission. We therefore consider two limiting values for this 
parameter f s = 0.5 and f s = 1. 

(v) We assume that the optical depth of the conical shell in the 
flat velocity range of the F-component is very high, such that it 
absorbs all the radiation of the central source it hides from our line 
of sight (in the flat part). 

(vi) We assume that the colliding winds shell has an hyperbolic 
shape. We consider its tilt due to the orbital motion. In calculating 
the shape and tilt the primary wind acceleration zone is considered. 

(vii) We assume, based on our previous papers (Kashi & Soker 
2008b,2009c), that the secondary is toward us near periastron (a> = 
90°). We take an inclination angle of i = 41°, and the other com- 
monly used binary parameters (semimajor axis, eccentricity, etc.). 

(viii) We note that after periastron passage some of the assump- 
tions break down, because accretion is likely to occur (Kashi & 
Soker 2009a; Akashi & Soker 2010). 

(ix) We calculate the part of the central source covered by the 
conical shell, for a) = 90° . 

(x) We calculate the covered part of the 1.083 /jm central con- 



tinuum emission source (equation [2}. By that we show that for 
to = 90° our model explains the observed increase in absorption 
strength of the F-component close to periastron passage (Fig. [5}, 
and explains its transient nature. 

(xi) We therefore conclude that the orientation is indeed that the 
secondary is toward us near periastron (a> = 90°). 

(xii) In addition, using hydrodynamic simulations (Fig.[2]and[3]> 
we find that the He I in the conical shell and in the pre-shocked 
secondary wind has a substantial column density in the velocity 
range -v = 1 150-1400 km s . We assume that the fraction (out of 
total He atoms and ions) of the He I in the 2 3 5 level is high enough 
to absorb most radiation in the -v < 1400 km s -1 (with a tail up to 
-v=* 1900 km s -1 ). 

The finding that the secondary is closer to us near periastron 
requires that some properties near periastron passage be explained 
by the accretion of the primary wind onto the secondary star (Kashi 
& Soker 2008a, 2009d). For accretion to occur onto the secondary 
star, that has a strong wind of his own, the binary must be close 
(Akashi & Soker 2010) and to interact strongly with the primary. 
The accretion that occurs near periastron passage is crucial to the 
understanding not only the present behavior of rj Car, but also its 
behavior during the major eruptions it undergone in the 19th cen- 
tury, the Great Eruption (GE) and Lesser Eruption. 

The debate on the absorbing source of the blue wing has im- 
plications far beyond the specific question on the orientation of the 
major axis of the binary system. The essence of the debate is the 
nature of the binary interaction process, which has implications on 
the nature of LBV major eruptions, and mass loss by very massive 
stars (e.g., Soker 2001, 2005, 2007; Kashi & Soker 2010; Smith et 
al. 2010; Smith 201 0a,b). 

The GE is the best studied example for a major LBV erup- 
tion, and serves as a test case for theories and models. A high rate 
accretion during the GE could have supplied the extra luminosity 
for 20 years, and the accreting secondary star could have launched 
two jets that shaped the bipolar structure - the Homunculus (Soker 
2001). Most likely other LBV major eruptions are also related to 
binary interaction (Kashi & Soker 2010), as was argued for the 
17th century eruption of P Cygni (Kashi 2010), rather than a sin- 
gle star phenomena as suggested by, e.g., Smith (2007). Smith & 
Owocki (2006) suggested that LBVs lose most of their envelope 
mass during major eruptions. Therefore, the presence of a strongly 
interacting companion can play a major role in the evolution of 
very massive stars. As suggested by Kashi et al. (2010), major mass 
transfer events in LBVs are related to optical transient objects and 
have a common powering mechanism - accretion onto a compan- 
ion star. Understanding the binary interaction in rj Car will shed 
light on other objects where binary interaction is thought to shape 
circumstellar nebulae, like planetary nebulae, symbiotic systems, 
and related objects such as the Red Rectangle. 
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APPENDICES (ASTRO-PH VERSION ONLY) 



APPENDIX A: THE BRIEF ABSORBTION IN THE 
CONICAL SHELL 

A support to a short duration absorption by the conical shell near 
periastron passage might come from the behavior of the Hydrogen 
HS line. The Hydrogen atom is a simple atom, lacking the com- 
plications resulting from absorption from a meta-stable state as for 
the He I /110830A line. The complication is in that it is not sim- 
ple to estimate the population of the relevant atomic state. Nielsen 
et al. (2009) present high spectral resolution (R ~ 80, 000) of the 
hydrogen Hc5 line, as well as of the He I /17067A line, across the 
2003.5 spectroscopic event. The He I /17067A line shows Doppler 
shift variations with the orbital phase that fit emission and absorp- 
tion in the acceleration zone of the secondary wind (Kashi & Soker 
2007, 2008b). The variations in the Doppler shift of the HS line 
are much smaller, perhaps ~ 0.3 times those of the He I /17067A 
line (hard to tell). This may suggest that part of the absorption in 
the HS line is from the conical shell, in addition to an absorption 
component associated with the primary wind. Close to periastron 
(phase 1.001 in figure 5 of Nielsen et al. 2009) the HS line shows 
an anomalous wide blue absorption wing, reaching velocities up to 
-v =i 900 km s -1 . At the closest observation times, phases 0.985 
and 1.115, this broad absorption does not exist. Though the blue 
wing is of slower velocities than observed in the He I /U0830A 
line, a resemblance between the two is very likely. Both lines show 
a high velocity absorption component just before periastron, which 
later disappears. Like the He I /U0830A line, it is possible that 
the H<5 absorption line originates in the conical shell, but in re- 
gions with higher density and lower velocity. We note that qual- 
itatively similar blue absorption wing close to periastron passage 
also appears in the Ha and Paschen rj lines, although less promi- 
nent (Davidson et al. 2005; for the Ha this trend is better seen in 
the animation prepared from these observations). Note that in the 
last periastron passage of 2009 the Ha (observed by Richardson et 
al. 2010) blue wing absorption component is not prominent, but it 
may exist and be difficult to see. 



APPENDIX B: ABSORPTION IN A JET? 

We here check whether an alternative absorbing material for the 
He I /110830A can be a collimated fast wind or a jet emitted from 
the secondary. Behar et al. (2007) analyzed CHANDRA MEG/HEG 
observations in the 3.7-7.5A band, where high excitation S, Si and 
Ar lines are observed. Far from periastron these lines show a sym- 
metric profile with no Doppler shift. Close to periastron the cen- 

troid moves to 200km s , and a wide blue wing is developed, 

reaching velocities of up to 2000km s -1 . Behar et al. (2007) in- 
terpreted the centroid motion to the orbital motion of the secondary, 
taking to = 90° . We note that though the observed velocity at phase 
-0.006 was out of their Doppler lines (their fig. 4), taking updated 
inclination angle of i = 41° instead of i = 53° as used by Behar 
et al. (2007) makes their fit much better and strengthen their claim 
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for the emission source for the centroids. The high velocity blue 
wing seen in these lines was attributed by Behar et al. (2007) to a 
collimated fast wind or a jet launched by the secondary. 

The wide blue wing at high velocities is clearly seen at phases 
-0.028(1 1.972) and -0.006(1 1.994) in fig. 3 of Behar et al. (2007), 
where the mean velocity profiles constructed from nine different 
spectral lines from five observations are plotted. The flux is stronger 
at phase -0.028 for velocities -v < 1400kms~' . Fig. 3 of Behar et 
al. (2007) presents normalized observation. From the original data 
presented in Henley et al. (2008) we learn that the true flux at phase 
-0.028 is ~ 1.7 times as strong as at phase -0.006 so that the wide 
blue wing at phase -0.028 is even more pronounced. Correcting 
for the true, not normalized, flux, we find that the wide blue wing 
at phase -0.028 is stronger than in phase -0.006 for not only for 
-v < 1400kms~' but up to -v < 2000kms~' . In this context, there 
is some difficulty in adopting the jet explanation for He I /U0830A 
line absorption at high blue velocities. At phase -0.028, 57 days 
before periastron, the accretion process have not yet started (Kashi 
& Soker 2009a), and hence an accretion disk around the secondary 
which is essential for producing jets does not exist. This is also 
supported by the X-ray flux which is strong and raising at phase 
-0.028, while weaker and decreasing at phase -0.006 (Corcoran 
2010), suggesting that substantial accretion has not yet started ac- 
cording to the accretion model (Soker 2005; Kashi & Soker 2009a). 
We therefore regard the explanation of a fast wind from a jet as the 
absorbing material for the He I /U0830A unlikely. 

It is possible that the gas emitting the highly ionized wide S, 
Ar, and and Si lines and the He I /U0830A absorption do have 
a common origin. The emission and absorption occur in different 
places along the outflow. It is therefore possible that the X-ray 
emitting gas is not a jet, but rather the secondary wind (pre and 
post-shock) itself. 



